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Ranking  and  Subset  Selection  Procedures  for  Exponential 
Populations  with  Type- I  and  Type- I I  Censored  Data 

by 

Roger  L.  Berger  and  Jee  Soo  Kin 

ABSTRACT 

* 

) 

In  this  paper  ranking  and  subset  selection  procedures  for 
exponential  populations  with  respect  to  the  largest  location  and 
scale  parameters  are  proposed.  The  data  are  assumed  to  be  gene¬ 
rated  from  Type- I  and  Type- II  censoring  mechanisms.  The  selec¬ 
tion  procedure  proposed  for  the  largest  scale  parameter  based  on 
Type- I I  censored  data  is  equivalent  to  Gupta's  procedure  for  gaoma 
populations.  Two  procedures  proposed  for  the  selection  of  the 
largest  location  parameter  under  Type- I  censoring  and  Type- I I 
censoring  are  asymptotically  equivalent. 
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1.  INTRODUCTION 

In  this  paper  we  investigate  the  problem  of  ranking  and 
selection  for  exponential  populations  with  Type- I  and  Type- I I 
censored  data.  We  consider  both  the  indifference  zone  formula¬ 
tion  of  Bechhofer  (1954)  and  the  subset  selection  formulation 
of  Gupta  (1956).  To  our  knowledge,  W.  Huang  and  K.  Huang  (1980) 
are  the  only  researchers  who  have  considered  a  selection  problem 
with  incomplete  data. 

Censored  or  incoaplete  data  arise  in  various  situations  such 
as  industrial  life-testing,  clinical  trials,  and  biological  expert 
ments.  To  motivate  our  discussion  and  to  illustrate  how  incom¬ 
plete  data  can  arise,  we  present  a  particular  situation  from  in¬ 
dustrial  life-testing.  A  batch  of  n  electronic  components  (or 
items)  are  placed  on  test  at  time  0,  and  the  experiment  terminates 
at  a  pre- specified  time  T.  The  failure  time  of  a  component  is 
obserwble  if  it  fails  before  time  T.  If  a  component  still  func¬ 
tions  at  the  closure  of  the  experiment,  its  failure  time  is  not 
observable.  The  component  is  said  to  be  censored  at  time  T. 

This  type  of  time  censoring  is  known  as  Type- I  censoring.  In 
this  life-testing  experiment  an  experimenter  might  not  know  be¬ 
forehand  what  value  of  the  fixed  censoring  time  T  is  appropriate, 
so  he  may  decide  to  continue  the  experiment  until  a  pre-specified 
number  r,  I  &  r  i  n,  (or  a  fraction  —,  or  a  proportion  100a%, 

0  <  a  <  1)  of  the  components  have  failed.  Those  components  still 
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functioning  at  the  tine  of  the  rth  failure  are  called  censored 
observations .  This  type  of  failure  censoring  is  known  as  Type- I I 
censoring.  Type- I  and  Type- I I  censoring  schemes  have  received 
nuch  attention  in  the  statistical  literature.  See  Bain  (1978), 

Barlow  and  Prcschan  (1966,  1967),  Epstein  and  Sobel  (1953),  and 
Mann  et  al.  (1974). 

We  deal  specifically  with  exponential  distributions  subject 
to  these  two  types  of  censoring  mechanisms,  location  parameter 
problems  are  considered  in  Section  2  under  Type-1  censoring  and  in 
Section  3  under  Type- I I  censoring.  In  Section  4,  ranking  and  se¬ 
lection  problems  for  scale  parameters  under  Type- II  censoring  are 
discussed.  Tables  for  implementing  the  new  procedures  are  presented. 
Throughout  this  paper,  we  use  the  word  increasing  in  place  of  non¬ 
decreasing  and  decreasing  in  place  of  nonincreasing. 
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2.  SUBSET  SELECTION  FOR  THE  LARGEST  LOCATION 
" ~*P ARSMBTeR'  1  bASHT  m  TVt>E-l~  @&5ftEiTbXW - 


Let  nx,  llj >  • ..»  nk  denote  k(k  2  2)  independent  exponential 
populations  with  density  functions 


f(x;  A.) 


V‘w*,/*  ,x4Xt 

0,  x  <  Aj,  i  «  1,  2,  ...»  k. 


The  scale  parameter  8  is  assumed  known.  The  location  parameter 
Aj  are  the  unknown  parameters  of  interest.  The  parameter  space 
will  be  denoted  by 


A  ■  (Aj ;  A^ »  (Aj,  A2<  ....  Xfc),  i  0). 

Let  A^j  £  £  ...  s  denote  the  ordered  values  of 

A}*  A2,  ...,  Xj,  and  let  denote  the  unknown  population  asso¬ 

ciated  with  X  ^ j .  Our  goal  is  the  subset  selection  goal  formu¬ 
lated  by  Gupta  (1956).  This  goal  is  to  select  a  nonenpty  sub¬ 
set  of  the  k  populations  containing  »  the  population  associated 
with  the  largest  parameter  A^j.  A  correct  selection,  denoted 
CS,  is  the  selection  of  any  subset  which  contains  The  popu¬ 

lation  is  referred  to  as  the  best  population.  If  more  than 
one  population  could  be  classified  as  best,  then  one  of  these 
is  arbitrarily  tagged  as  the  best.  This  is  done  only  for  the 
purpose  of  evaluating  the  infimum  of  the  probability  of  a  CS. 

We  shall  propose  a  selection  rule  which  satisfies  tho  P*-condition, 


I 

! 
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inf  P.  (CS|R.)  fc  P*,  (2.1) 

A  —  1 

where  -■  <  P*  <  l  is  specified  by  the  experimenter  prior  to 
experimentation . 

For  i  •  1,  2,  k  let  X^,  j  ■  1,  2,  ....  n  denote  a 
random  sample  of  size  n  from  n^.  The  X^.  might  be  the  life 
lengths  of  n  items  from  population  placed  on  test  at  time  0. 

The  X_  are  not  themselves  observable.  Rather  they  are  subjected 
to  Type-1  censoring.  That  is,  there  is  a  specified  time  T  >  0 
past  which  no  observation  will  take  place.  Thus  the  observable 

data  are  i  »  1,  2 . k,  j  •  1,  2,  ....  n,  where  X*^  ■  minfXjj,  T). 

Using  the  Factorization  Theorem  (Lehmann  (1959) ,  P.  49)  it  can 
be  shown  that  (Yj,  Y2>  ....  Y^)  is  a  sufficient  statistic  in 
this  problem  where 


Ya  »  min(X*lt  ....  X*n)  -  min(min(X.  j. 


Xin>’ 


T). 


The  selection  procedure  we  propose  depends  on  the  data  only  through 
the  sufficient  statistic  (Yj,  Y2,  ...»  Y^). 

The  observations  Yj,  Y2,  . ...  Y^  are  independent  random 
variables.  The  c.d.f.  of  Y^  is  G(y^;  Ai)  where 


G(y;  X) 


0 

1  .  e-n(y-X)/e 
1 


if  y  <  X 
if  X  s  y  <  T 
if  T  sy 


if  X  <  T  and 


G(y;  X) 


0  if  y  <  T 
1  if  T  sy 


1 


4 
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if  T  s  X.  If  X  <  T,  G(y;  X)  is  a  mixed  probability  distribution 
with  an  absolutely  continuous  part  defined  by  the  density 
g(y;  X)  »  £  e'n(y-X)/6  £or  x  i  y  <  T  and  a  discrete  part  con¬ 
sisting  of  a  point  mass  at  T  with  probability  e“n(T"*)/®.  If 
T  s  X,  G(y;  X)  is  degenerate  at  T. 

The  selection  procedure  we  propose  is  denoted  by  Rj. 

Procedure  Select  JI^  if  and  only  if 


Y4  *  max  Y.  -  22. 
1  lijSk  *  n 


(2.2) 


where  d  *  d(k,  P*.  0,  T,  n)  4  0  is  chosen  so  the  P*-condition 
(2.1)  is  satisfied. 


The  method  of  computing  d  will  be  given.  It  will  be  shown 

that  d  depends  on  n,  T  and  e  only  through  —  and  that  0  s  d  5  ^  . 

nT 

Table  1  gives  values  of  d  for  k  *  2(1)S,  —  ■  ,25(.25)4.00  and 

5(1)10  and  P*  «  .75(.05).95  and  .99. 

If  d  a  $  then  max  Y.  -  ^  s  0  with  probability  1  for  any  X. 

8  lsjsk  3  n 

Thus  Pa(CS|Rj)  »  l  for  all  V.  For  P*  <  1,  (2.1)  can  be  achieved 
with  a  d  <  -jp. 

We  will  now  examine  Px(CS|Rj)  to  determine  how  d  should  be 

chosen  to  ensure  the  P*-condition  (2.1). 

.  d0 

Theorem  2.1.  Let  d  ■  — .  Then 

1  for  T  -  J  s  lw 

/*'3  B.  e-n()"*[k])/9dy  (2.3) 

*P0  3*1 

♦  e"n(T"J"V],/9  for  i[k]  <  T  -  d. 
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Proof.  Let  i  *  1»  2,  ....  k  be  the  observation  from 


or 


if  T  -  I  SX[k],  i  -  Px(Ym  a  T  -  d)  S  Pxtt^  *  max  Ym-H) 


-  PxCCSlRj)  «  l. 


‘lv‘00 


1-00  '  ^  *0) 


If  *w  <  T  -  d. 


PxlCSlRj)  *  tV(Y(k)  *  Y(j)“d.  j-1,  2,  ....  k-1,  and  Y  <  T  -  d) 

4  *VCY(k)  *  Y(j)-^’  3*1'  2*  •**»  k"1»  md  Y(k)  * T  ■  d> 

«  Px(Y(lc)  4  d  *  Y(j)»  3"1*  2‘  •***  k'l»  mi  Y(k)  <  T  -  d) 


4  P^oo  *  T  *  d> 

rT-J  k-1 

»  n  G(y  ♦  d;  Xm)g(y;  Xrkl)dy 

m  !kl 

♦  o-Cr-^x[k])/8 


Substituting  for  G  and  g  yields  (2.3).  || 

For  a  fixed  value  of  X^ ,  the  P^  (CS | R^)  in  (2.3)  decreases 
if  k[2]*  •••»  ^^k]  are  a11  replaced  by  X^j.  Thus  the 

inf  Px^lRj)  takes  place  when  Xjjj  ■  Xj2j  *  ... 

X  •  (X,  X,  ...,X)  and  X  <  T  -  1,  then 


x[k)*  If 


Px^lRj)  ■  J*  d{l  -  e"n(y4d’X)/tl)k’1  2.  e-n(y-X)/8  dy 


*  e-n(T-J-x)/e 

-  IedU  -  e-"(y^-X)/8)k|Y-a'  „  e-n(T-I-X)/0 

-  i  edU  -  e-n(T‘X)/9}k  -  £  ed(l  -  e^  ♦ 
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Differentiation  with  respect  to  X  shows  that  P^(CS|Rj)  is  a 
nonotone  increasing  function  of  X  for  0  *  X  <  T  -  d . 

Therefore,  P^(CSlRj)  is  minimized  by  setting  X  o  o. 

Theorem  2.2.  For  any  value  of  d,  0  £  d  *  — , 

.»T  .51 

Inf  P^CCSlRj)  •  ed(i(l  -  e  "*")*  -  J-(1  -  e'd)k  .  e  9  >.  (2.4) 


We  have  restricted  d  to  be  nonnegative  so  that  a  nonempty 


subset  is  always  selected.  If  d  ■  0,  then  from  (2.4)  we  see  that 

nT  -nT 

inf  P*  (CS|  R.)  ■  r<l  -  e)  6  )k  *  e  If  P*  is  less  than  this 
A  -  1  “ 

value,  then  d  ■  0  can  be  used  and  strict  inequality  holds  in  (2.1). 
This  can  be  observed  for  some  small  P*  values  in  Table  1.  But  if 
P*  is  greater  than  or  equal  to  this  value,  d  can  be  chosen  so  that 
equality  holds  in  (2.1)  by  equating  the  right  side  of  (2.4)  to  P* 
and  solving  for  d.  Note  that  (2.4)  depends  on  n,  6,  and  T  only 


nT 


through  -g-.  Thus  we  have  Theorem  2.3.  nT 


nT 
'  9 


Theorem  2.3.  For  k  i  2,  ~  >  0  and  J(1  -  e  6  )k<  e  s  s  P*  <  1, 
if  d  *  d(k,  P#,  »-ln  u  where  u  is  the  solution  of  the  polynomial 


nT 


nT 


(1  -  u)k  ♦  kP*u  -  ke  0  -  (1  -  e  0  )k  »  0  (2.5) 

.nT 

in  the  interval  e  <  u  £  1,  then  the  P*-condition  (2.1)  is  satis¬ 
fied  with  equality. 

The  values  of  d  in  Table  1  were  confuted  by  solving  (2.5). 

Let  0  »  (0,  0,  ....  0).  Since  P0(CS|Rj) (the  expression  in  (2.4)) 


a 
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is  a  continuous  strictly  increasing  function  of  d  for  0  <  d  s  ~ 
and  P0(CS |Rj)  ■  1  if  d  ■  (2.5)  always  has  a  unique  solution 

nT 

in  e  9  <  u  s  1.  As  ~  +  <*,  d(k,  P* ,  —)  d(k,  P#,  *)  where 
d(k,  P *,  •)  is  the  solution  of 

ed{i  -  £<1  -  e‘d)k>  «  P*.  d  2  0.  (2.6) 

Equivalently  d(k,  P*,  »)  a  -in  u  where  u  is  the  solution  of 

(1  -  u)k  ♦  kP*u  -  1  »  0,  0  <  u  <  1. 

The  values  of  d(k,  P*,  ®)  are  also  listed  in  Table  1.  For  the 
values  of  k  and  P*  in  Table  1, 

|d(k,  P*.  •)  -  d(k,  P*,  10)|  s  .0001. 

Finally,  we  note  sane  reasonable  properties  of  Rj  in  Theorem 

2.4. 

Theorem  2.4.  Fix  T  >  0  and  e  >  0. 

(i)  If  T  s  or  <  Ajkj  <  T.  lim  P\  CCSlRi)  •  1. 

(ii)  (Monotonicity)  If  Ai  s  Aj,  (select  IIjlRj)  s  Px(select  jijlRj)* 

(iii)  If  <  «»inimum(Ajkj,  T),  lim  Ex(S|Rj)  ■  1 

where  S  is  the  site  of  the  selected  subset. 
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Proofs 

(i)  If  T  s  A^,  ^(CSjRj)  ■  1  *OT  every  n.  Suppose 

l(k-l]  *  X[k]  *  T‘  Since  r(t)  ■*  X(i]  *■*• 

PjCCSl^)  2  ^  Y(J))  •  1  u  n  *  •< 

(ii)  G(y;  X)  is  stochastically  increasing  in  x,  i.e.. 


G(y;  Aj)  s  G(y;  A^  if  A^  s  A^ .  This  fact  and  Problem  11,  page 
112  of  Lehmann  (1959)  yield 


/T  k 

Px(select  HjRj)  «  n  C(y  ♦  J  j  Affl)dG(y;  V 
—  *  0  m»l  1 


i 


•T  k 


de 


5  n  G(y  ♦  ^  ;  A  )dG(y;  A  ) 
’0  m«l  n  “  x 

m*j 


fT  k  de 

H  G(y  ♦  ;  A )dG(y;  A  ) 

*0  m*l  n  J 


m*j 


■  P^Cselect  n^lRj). 


(iii)  Ex (S|Rl) 


k 

■  Jj  V*lw 


By  (i)  and  (ii)  it  suffices  to  show  lim  P^ (select  n0c_i)lRi)  *  0 


Since  lim  d  -  lim  d(k,  P*.  ~)  ■  d(k,  P*f  •)  <  •,  lim  -  0. 


de 


in 


n-H» 


n 


-  f-  < 
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3.  SUBSET  SELECTION  FOR  THE  LARGEST  LOCATION 
PARAMETER  BASED  ON  TYPE- I I  CENSORED  DATA 


Assume  that  for  i  »  1.  2,  ....  k,  X^,  )  »  1,  2,  ...»  n  de¬ 
note  a  random  sample  of  size  n  from  n. .  Thus  the  X. .  are  defined 

ii) 

as  in  Section  2.  Again,  we  are  interested  in  selecting  a  subset 

containing  A^j.  In  this  section,  we  consider  Type-II  censored 

data,  rather  than  the  Type- I  censored  data  of  Section  2. 

Let  r^,  i  -  1,  ....  k  be  fixed  integers  satisfying  1  s  ^  i  n. 

Let  X.  ri,  £  X, s  ...  s  X.  r  ,  denote  the  order  statistics  from 
X[1J  i 12J  ilnj 

population  n. .  The  data  from  population  is  Type- I I  censored 

if  the  observations  available  from  n4  are  X. ril,  . ..,  X. .  , .  This 

i  1UJ  1  lxi^ 

type  of  data  arise  if  n  items  are  placed  on  test  simultaneously 
and  observation  stops  afteT  the  i\ tb  failure. 

Assume  that  Type- 1 1  censored  data  are  available  from  each 
IT,  i  *  1,  2,  ...,  k,  i .  c . ,  the  data  are  X^  ^  j  j ,  z  *  1,  2,  •  •  * ,  k. 


j  •  1,  2,  ij.  Using  the  Factorization  Theorem  (Lehmann,  1959), 
it  can  be  shown  that  (X^^,  X2[i]’  *•*»  ^Scfi]^  is  a  sufficient 

statistic  for  this  problem.  Thus  the  experimenter  needs  only  ob¬ 
serve  the  first  failure  time  from  each  population.  Alternatively, 


(X1  (11 '  X2[l]’  *“*  ^[l]5  “X  be  thou8ht  of  83  the  result  of 

Type-II  censoring  with  r^  »  1,  i  ■  1,  ...»  k. 

Ofosu  (1974)  proposed  the  following  subset  selection  rule  for 

x(k]  b8Sed  0,1  (xi[iK  x2(i]’  Xkfl]) ' 
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Procedure  R^:  Select  if  end  only  if 

Xltn  2  "idi  . k 

where  d  »  d(k,  P*)  i  0  satisfies 

£  ed(l  -  (1  -  e"d)k)  -  P*  (3 

Ofosu  showed  that  for  £  <  P*  <  1,  satisfies  the  P*-condition 

inf  PxCCSjRj)  i  P* 

if  d  is  chosen  as  in  (3.1).  Condition  (3.1)  is  the  limiting 
case  (2.6)  for  the  procedure  Rj  discussed  in  Section  2.  This  is 
unsurprising  since  1^0^  -  X.  ^j)  ■*  1  as  T  for  fixed  n  and  0. 

Some  desirable  properties  of  the  rule  R2  are  listed  below. 
Properties  (i)  and  (ii)  were  noted  by  Ofosu  (1974) .  Properties 
(iii)  and  (iv)  can  be  proven  in  a  similar  fashion  to  Theorem  2.4. 

(i)  sup  E.  (S|R,)  •  kP\ 

(ii)  If  Xj  s  P^ (select  n.jR2)  <  P^ (select  n^l^)* 

(iii)  If  <  A,kj,  lim  Px(CS|R2)  *  1. 

(iv)  If  *  X  [kj »  11®  ^ 

Property  (i)  is  in  marked  contrast  to  Rj  in  the  Type- I  cen¬ 

soring  model.  If  2  T,  Ex(S|Rj)  ■  k  (all  the  populations  are 
surely  selected).  Berger  (1979)  showed  that  if  a  selection  rule 


satisfies  (i)  then  the  rule  is  ainiaax  with  respect  to  S  and  S', 


the  number  of  non-best  populations  selected.  That  is,  satisfies 


sup  E, (S Jr.)  ■  inf  sup  E, (S|R) 
A  ^  2  V  .  A  i 

P 

and 


sup  E.(S'|R,)  -  inf  sup  E.(S'|R) 

A  *  2  V  *  A  - 

P 

where  D  *  is  the  set  of  all  rules  which  satisfy  the  P*-condition. 
P 

In  the  Type- I  censoring  model  of  Section  2,  no  nonrandomized  rule 
can  be  minimax  since  P^O^  *  T,  i  *  1,  2,  ...»  k)  •  1  if  T. 
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4.  RANKING  AND  SUBSET  SELECTION  FOR  SCALE  PARAMETERS 
BASfeb  ON  TYPE- IltfeNSORED  DAtA. 


Let  IIj,  n2,  . ...  be  k(k  2  2)  exponential  populations  with 
density  functions 


f(x;  6^  . 


^  1^  e’x^®i  *  **  x  *  0  ^  «i  >  0 


(4.1) 


otherwise. 


i  *  1,  2,  ...,  k. 


The  ordered  values  of  the  unknown  scale  parameter  are  denoted  by 


9[1]  5  e{2]  *  “*  S  0[k]‘  (4,2) 

The  unknown  population  associated  with  8^  is  denoted  by  • 

It  is  assumed  that  there  is  no  prior  knowledge  about  the  correct 
pairing  of  the  k  populations  and  the  ordered  scale  parameters  8^. 
In  the  context  of  the  reliability  theory  and  life-testing  models 
is  the  mean  life  of  an  item  whose  life  length  is  described  by 
an  exponential  distribution  function.  The  population  with  the 
parameter  equal  to  6^  is  called  the  best  population.  The  goal 
is  to  identify  the  best  population,  and  also  the  t(l  s  t  a  k  -  1) 
best  populations. 

From  each  population  n^,  i  ■  1,  2,  ...,  k  we  take  a  sample 

of  n  items.  Let  s  Xi[2]  s  s  Xi[n]  <*enote  the  or<Ier  sta¬ 

tistics  representing  the  failure  times  of  the  n  items  from  popula¬ 
tion  ni.  We  consider  Type- I I  censoring  as  in  Section  3.  Let  r 
be  a  fixed  integer  such  that  1  s  r  s  n.  Under  Type- I I  censoring 


the  first  r  failures  from  each  population  are  to  be  observed. 
Observation  on  ceases  with  the  observation  of 
(n  -  r)  items  whose  failure  times  are  not  observable  be¬ 
yond  the  Xi[rJ  become  censored  observations.  Thus,  the  length 
of  the  time  population  IT  is  observed  is  the  random  variable 
Xi[r]*  1VPe“11  censoring  was  first  investigated  by  Epstein  and 
Sobel  (1953). 

Section  4.1  introduces  a  statistic  known  as  the  total  time 
on  test  (TTOT)  statistic  generated  from  Type- I I  censoring.  Some 
of  its  well-known  properties  are  stated  in  Section  4.1.  Selec¬ 
tion  rules  are  proposed  on  the  basis  of  the  TTOT  statistics. 

4.1.  Total  Time  on  Test  (TTOT)  Statistics  and  Its  Properties. 

For  the  Type- II  censored  sampling  the  likelihood  function 
of  the  first  r  failures,  X^j,  X^j,  ••••  Xi[r]  from  n*  is 
given  by 


f (xi[l] ,Xi[2] • *  *  * »xi[r] ; V 


,  .  r 

- r  exp[-¥  (  l  x.f. i*(n-r) x.. 

(n-r)/0[  ei  jil  iI3J  il 


for  0  s  s  xiJ2J  s  ...  sxiIrj 


0,  otherwise. 


(4.: 


The  maximum  likelihood  estimate  (MLE)  0^  of  0^  obtained  by  setting 
d  Iok  f  _ 


A  1  _ 
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* 

where  •  l  xjjjj  ♦  (n  -  r)x^jrj*  which  can  be  written  as 

Zi  *jH"  -  i  •  *>^01  -  X1[M]>  ♦  <"  -  r>xi[r]'  {4  S> 


where  X^qj  =  0.  represents  the  total  tine  spent  on  observing 

n  items  undergoing  life-testing  until  the  rth  failure  from 

whether  an  item  is  censored  or  uncensored.  Z^  is  called  the  TTOT 

* 

statistic.  We  state  some  of  the  known  properties  of  See, 
for  example,  Epstein  and  Sobel  (1953),  Mann  et  at.  (1974)  and 
Bain  (1978). 

Lemma  4.4.1.  defined  in  (4.4)  has  the  following  properties. 

A 

(i)  8^  is  a  complete  and  sufficient  statistic  for  e^. 


(ii)  is  the  MLE  of  8^ 


(iii)  6.  is  the  UMVUE  of  8^ 

(iv)  8^  is  a  strongly  consistent  estimator  of  e^,  if  —  is 
fixed. 

2r®i  2 

(v)  -s —  is  distributed  as  x  ,a  chi-square  with  2r  de- 
9i  (2r) 

grees  of  freedom. 

A 

(vi)  Zi  •  r8^  has  a  gamma  distribution  with  shape  parameter 
r  and  scale  parameter  8^. 

Remark  4.1.1.  The  distribution  of  Z^  is  independent  of  the  sample 
size  n.  The  advantage  of  r  failures  out  of  a  sample  of  n  items 
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(n  >  r)  over  r  failures  out  of  a  sauq^le  of  r  items  is  in  the 

savings  of  time.  The  expected  waiting  time  to  observe  the  rth 

r 

failure  in  the  sample  of  n,  E(Xijrj)  *  I  JpJ+T  *  is  less 

than  the  expected  waiting  time  for  all  r  failures  in  the  sample 
of  r  items. 

A 

Procedures  based  on  Z^,  equivalently  on  8^,  are  proposed 
with  both  formulations,  indifference  tone  formulation  and  sub¬ 
set  selection  formulation,  in  mind. 

4.2.  Indifference  Zone  Formulation. 

Our  goal  is  to  select  t,  1  £  t  s  k  -  1,  populations  asso¬ 
ciated  with  the  t  largest  scale  parameters  e^j,  8  »  k  -  t  ♦  1, 
k  -  t  ♦  2,  ....  k  whenever  the  unknown  parameter  configuration 
j8  »  (0j,  e2,  . ...  0^)  is  assumed  to  lie  outside  some  subset  of 
the  entire  parameter  space,  denoted  by  e.  The  selection  of  the 
subset  containing  exactly  the  t  largest  populations  will  be  called 
the  correct  selection  (CS).  The  preference  zone,  0(6)  is  defined 
by 

0(6)  «  {£  e  0:  6[k-t]/e[k-t*l]  <  «>  (4.6) 

where  0  <  6  <  1  is  fixed  by  the  experimenter. 

We  shall  propose  a  selection  procedure  Rj  for  which  the  pro¬ 
bability  of  the  correct  selection  satisfies  the  following  inequality 

P  (CSlR.)  iP*.  -~  <  P*  <  1  for  all  6  c  0(6). 

-  3 


(4.7) 
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where  6  and  P*  are  specified  in  advance.  This  is  the  formula¬ 
tion  of  the  selection  problem  presented  by  Bechhofer  (19S4). 

Let  the  ordered  values  of  be  denoted  by 

Z{1]  5  Z[2]  S  *•*  *  Z[k]  (4*8) 

The  procedure  Rj  is  defined  as  follows. 

Procedure  R^:  Select  t  populations  corresponding  to  Z^-j , 

Z[k-1]’  *"*  Z[k-t*l]  ^i‘e*  select  P<>P«la-  (4.9) 
tion  for  which  Z^  2: 

The  rule  R^  which  depends  on  observations  (censored  or  uncensored) 
only  through  the  TTOT  statistics  Z^  is  equivalent  to  the  procedure 

A 

based  on  the  MLE's  8^. 

Having  specified  two  quantities  6  and  P*  an  experimenter  is 
now  interested  in  determining  the  ams2iest  value  of  r  to  guarantee 
the  basic  probability  requirement  (4.7).  For  this  purpose  the 
infimum  of  (CS| R^)  must  be  obtained.  Let  the  TTOT  statistic 
from  the  population  associated  with  8^j  be  denoted  by  Z^j .  It 
can  be  shown  that 

P0 (CS|  Rj)*Pq (®ax(Z (1j , . . .  »Z(]{_t^  <  ®^n  ^Z(k-t*l)  *  * "  ’  *Z(k)  ^  (4.10) 

and  that  the  least  favorable  configuration  in  6(6)  is  given  by 

®[1]  *  ®[2]  *  •**  "  ®(k-t]  *  50(k-t*l]’  and  ®[k-t*l]"  *  9[k)- 

Thus,  we  can  obtain 

inf  Pfl(CS|R  )  •  (k-t)|F*'t'1(u)(l  -  F  (gu))tdF  (u)  (4.11) 

6(6)  -  5  *0  T  r  r  » 


x  e  dx. 


where  Fr(u)  -  m  Jq 

Since  2 ^ are  independent  gamma  variables  with 
common  known  shape  parameter  r  and  unknown  scale  parameter  0^, 
the  selection  rule  is  closely  related  to  the  selection  rule 
studied  by  Bechhofer  and  Sobel  (1954)  in  their  investigation  of 
the  ranking  variances  of  normal  populations  based  on  the  sample 
variances.  Tables  in  Bechhofer  and  Sobel  (1954)  or  Gibbons,  Olkin 
and  Sobel  (1977)  can  be  used  to  determine  the  minimum  r  required 
to  satisfy  (4.7). 

We  now  study  the  behaviour  of  the  P0(CSfR3)  as  6  and  n  vary. 
Theorem  4.2.1.  For  any  fixed  r,  1  s  r  £  n, 

lim  inf  P  (CS | R_)  «  1. 

<5+0  0(6)  -  3 

Proof.  This  follows  readily  from  the  expression  (4.11), 

inf  Pe(CS|R  )  -  (k  -  t)  f  F*'t_1(u)  (1  -  FT(6u))tdF  (u) 

0(6)  -  3  >o  T  r  r 

♦  (k  -  t)  I  F^_tl (u)dF(u)  as  6  +  0 
>0  T  r 

•  I-  II 

Suppose  that  the  experimenter  specifies  a  fraction  p  say 
q(0  <  q  <  1),  instead  of  fixing  r.  Then,  as  a  consequence  of 
Lemma  4.1.1  (iv)  we  get: 
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Theorem  4.2.2.  Let  ^  «  q  be  a  fixed  fraction.  Then  for  any 
8  c  0(6).  lim  P0(CS|R3)  «  1. 


4.3  Subset  Selection  Formulation 

The  approach  we  take  in  the  present  section  is  to  select  a 
subset  of  k  given  populations  which  contain  the  population  asso¬ 
ciated  with  the  largest  scale  parameter  A  correct  selection 

(CS)  is  defined  as  the  selection  of  any  subset  which  contains  the 
best  population.  We  require  that  the  probability  of  including 
the  best  population  in  the  subset  selection  to  be  at  least  a 
p re-assigned  fixed  number  P*,  that  is. 


inf  P  (CS | R)  2  P*,  where  f  <  P*  <  1. 

a  9  K 


e  * 


(4.12) 


4.3.1  Selection  Procedure  R4  and  Its  Properties. 

Based  on  the  TTOT  statistics  2.  *  f  X, ...  ♦  (n  -  r)X 

i  »  1,  2,  . ...  k,  W.  Huang  and  K.  Huang  (1980)  proposed  the  follow¬ 


i[r]  ’ 


ing  procedure. 

Procedure  R^:  Select  the  population  if  and  only  if 


Zi  *  c  max  Z.  (4.13) 

1  lSjSk  3 

where  c  is  a  positive  constant,  0  <  c  <  1,  pre-assign ed  so  as  to 
satisfy  the  P*-condition  (4.12). 

Since  Z^  has  a  gamma  distribution  with  the  known  shape  para¬ 
meter  r  and  the  scale  parameter  8^  (Lemma  4.1.1  (vi)),  the  expres¬ 
sion  for  the  inf  PA  (CS|  R.)  and  all  the  desirable  properties  of  the 
e  _  4 


■J 


rule  R4  are  precisely  the  same  as  those  of  the  rule  proposed 
by  Gupta  (1963).  Gupta  (1963)  considered  a  problem  of  a  selec¬ 
tion  and  ranking  procedure  for  k  gamma  populations  with  respect 
to  the  largest  scale  parameter  when  the  shape  parameter  are 
assumed  to  be  known  and  all  equal.  He  derived 

inf  P0(CSJR4)  -  jV'^dF^u)  (4.14) 

where  Fy(u)  -  ~  J  xr'Vxdx. 

We  note  that  the  right-hand  side  of  (4.14)  is  independent  of  the 

sample  size  n.  It  depends  only  on  r,  the  number  of  observed  failure 

times.  For  various  values  of  k,  r  and  P*,  the  associated  constants 

c  satisfying  the  P*-condition  are  tabulated  in  Gupta  (1963). 

As  a  possible  question  in  designing  experiments  one  may  ask: 

For  a  fixed  value  of  the  constant  c,  0  <  c  <  1,  what  should  the 

amount  of  the  observed  failures  in  the  data  be  so  as  to  satisfy 

the  pre-assigned  probability  P*  specified  by  the  experimenter? 

Monotonicity  of  the  inf  P  (CS|R.)  will  resolve  this  problem  and 

9  1 

will  give  us  a  unique  smallest  integer  r  for  which  inf  P. (CS |RJ  *  P  . 

.  ©  1 

Theorem  4.3.1.  For  any  given  values  of  k,  c  and  P  in  the  selec¬ 
tion  rule  R4 

inf  Pe(CS|R4)  -  £  1 (~)dFy (u) 

where  Fr(u)  •  j^tr“1e“tdt.  is  strictly  increasing  in  r  >  l. 
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Proof.  Suppose  that  X  is  a  random  variable  having  the  density 
function  gy(x)  ■  cr  exp[(l  -  c)F^(x)]  for  x  c  {0,  1]  and  r  >  1 
where  F~*(x)  is  the  inverse  function  of  the  incomplete  gamma, 

Fr(u)  ■  |^yT‘Vydy.  Put  B(x)  -  log  gr^(x)  -  log  gr^(x) 

for  r2  >  Tj  >  1  and  0  i  x  s  1. 

— a - a-c).-.A(xl - -  where  f  (x)  is  a  gamma  density 

*r  [Fll(x)]f  [F"1 (x) ]  r 

T1  T1  r2  *2 

and  A(x)  *  f  [F^Cx)]  -  f  (F-1(x)].  Alam  (1970)  shows  that  A(x) 
rl  rl  r2  r2 

is  nonnegative  for  all  x,  x  €  [0,  1].  This  proves  B(x)  is  in¬ 
creasing  in  x  for  all  c,  which  in  turn  implies  gr(x)  is  a  density 

function  having  a  monotone  likelihood  ratio  in  x.  Hence, 
k- 1  f *  k- 1 

E  (X  )  ■  x  g  (x)dx  is  increasing  in  r  by  Lehmann's  lemma 
*0  r 

(19S9,  page  74).  But 

Er(Xk_1)  -  £  F^’1(H)dFr(u). H 

Suppose  that  the  experimenter  observes  a  fraction  £  of  failures 
Then  the  procedure  is  a  consistent  procedure  in  a  sense  that 
for  a  sufficiently  large  n,  the  probability  of  including  the  best 
population  in  the  subset  selected  approaches  1. 

Theorem  4.3.2.  Let  £  ■  q,  0  <  q  <  1,  be  fixed.  If  ®[k-l]  *  ®[k]' 

then  lim  P0(CS|R  )  ■  1. 
ir*«  -  4 

Proof.  This  is  immediate  from  Lemma  4.1.1  (iv),  because  R4  is 
based  on  Z. ,  or  equivalently  on  9.  defined  in  (4.4).  || 
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Table  1.  This  table  gives  the  necessary  d- value  required 
for  the  procedure  R^. 


2 

P*  a 

3 

.75 

4 

S 

-5 

.25 

.OOOO1 

.OOOO2 

.OOOO3 

.oooo4 5 

.0000 

.50 

.0987 

.1818 

.2030 

.2093 

.2123 

.75 

.2428 

.3847 

.4297 

.4477 

.4623 

1.00 

.3642 

.5681 

.6405 

.6736 

.7123 

1.25 

.4616 

.7309 

.8342 

.8851 

.9623 

1.50 

.5360 

.8711 

1.0090 

1.0805 

1.2123 

1.75 

.5898 

.9872 

1.1628 

1.2577 

1.4623 

2.00 

.6268 

1.0789 

1.2934 

1.4144 

1.7123 

2.25 

.6514 

1.1478 

1.3996 

1.5482 

1.9623 

2.50 

.6672 

1.1970 

1.4820 

1.6578 

2.2123 

2.75 

.6772 

1 . 2305 

1.5426 

1.7433 

2.4623 

3.00 

.6834 

1.2526 

1.5851 

1.8068 

2.7123 

3.25 

.6872 

1.2669 

1.6137 

1.8516 

2.9623 

3.50 

.6895 

1.2758 

1.6324 

1.8820 

3.2123 

3.75 

.6909 

1.2813 

1.6442 

1.9018 

3.4623 

4.00 

.6918 

1.2847 

1.6517 

1.9145 

3.7123 

5.00 

.6930 

1.2893 

1.6619 

1.9324 

4.7123 

6.00 

.6931 

1 . 2900 

1.6633 

1.9349 

5.7123 

7.00 

.6931 

1.2900 

1.6635 

1.9352 

6.7123 

8.00 

.6931 

1.2901 

1.6636 

1.9353 

7.7123 

9.00 

.6931 

1.2901 

1.6636 

1.9353 

8.7123 

10.00 

.6931 

1.2901 

1.6636 

1.9153 

9.7123 

00 

.6931 

1.2901 

1.6636 

1.9353 

00 

1 Actual  inf  P(CS| Rj)  ■  .8033 

2 Actual  inf  P(CS| Rj)  ■  .7824 

3Actualinf  PCCSlRj)  »  .7794 

4 Actual  inf  P(CS|Rj)  •  .7789 

5d  ♦  ~  ♦  In  P*  as  k  -►  - 
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P* 

o 

00 

« 

nT\  . 

T\ 

2 

3 

4 

5 

OD 

.  25 

.OOOO1 

.0222 

.0261 

.0267 

.0269 

.50 

.1760 

.2497 

.2683 

.2740 

.2769 

.75 

.3383 

.4610 

.4988 

.5141 

.5269 

1.00 

.4801 

.6566 

.7168 

.7442 

.7769 

1.25 

.5989 

.8350 

.9210 

.9629 

1.0269 

1.50 

.6939 

.9938 

1.1098 

1.1687 

1.2769 

1.75 

.7658 

1.130S 

1.2806 

1.3595 

1.5269 

2.00 

.8176 

1.2435 

1.4309 

1.5329 

1.7769 

2.25 

.8531 

1.3325 

1.5584 

1.6862 

2.0269 

2.50 

.8766 

1.3989 

1.6617 

1.8169 

2.2769 

2.75 

.8917 

1.4462 

1.7415 

1.9237 

2.5269 

3.00 

.9011 

1.4783 

1.8000 

2.0067 

2.7769 

3.25 

.9070 

1.4994 

1.8408 

2.0680 

3.0269 

3.50 

.9106 

1.5129 

1.8682 

2.1111 

3.2769 

3.75 

.9129 

1.5214 

1.8860 

2.1402 

3.5269 

4.00 

.9142 

1.5266 

1.8973 

2.1591 

3.7769 

5.00 

.9160 

1.5338 

1.9130 

2.1864 

4.7769 

6.00 

.9163 

1.5347 

1.9152 

2.1902 

5.7769 

7.00 

.9163 

1.5349 

1.9155 

2.1908 

6.7769 

8.00 

.9163 

1.5349 

1.9156 

2.1909 

7.7769 

9.00 

.9163 

1.5349 

1.9156 

2.1909 

8.7769 

10.00 

.9163 

1.5349 

1.9156 

2.1909 

9.7769 

OB 

.9163 

1.5349 

1.9156 

2.1909 

OO 

1 Actual  inf  P(CS|Rj)  ®  .8033 
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P*  -  .85 


f\‘ 

2 

3 

4 

5 

- 

.25 

.0586 

.0831 

.0867 

.0873 

.0875 

.50 

.2551 

.3151 

.3302 

.3350 

.3375 

.75 

.4376 

.5355 

.5651 

.5771 

.5875 

1.00 

.6028 

.7441 

.7907 

.8118 

.8375 

1.25 

.7478 

.9394 

1.0060 

1.0381 

1.0875 

1.50 

.8702 

1.1189 

1.2095 

1.2547 

1.3375 

1.75 

.9687 

1.2799 

1.3988 

1.4598 

1.5875 

2.00 

1.0440 

1.4196 

1.5714 

1.6512 

1.8375 

2.2S 

1.0986 

1.5360 

1.7243 

1.8263 

2.0875 

2.50 

1.1363 

1.6285 

1.8551 

1.9822 

2.3375 

2.75 

1.1614 

1.6982 

1.9621 

2.1162 

2.5875 

3.00 

1.1775 

1.7481 

2.0454 

2.2267 

2.8375 

3.25 

1.1877 

1.7823 

2.1070 

2.3133 

3.0875 

3.50 

1.1940 

1.8048 

2.1503 

2.3778 

3.3375 

3.7S 

1.1979 

1.8192 

2.1795 

2.4236 

3.5875 

4.00 

1.2003 

1.8283 

2.1985 

2.4545 

3.8375 

5.00 

1.2035 

1.8409 

2.2258 

2.5010 

4.8375 

6.00 

1.2039 

1.8426 

2.2297 

2.5079 

5.8375 

7.00 

1.2040 

1.8429 

2.2302 

2.5088 

6.8375 

8.00 

1.2040 

1.8429 

2.2303 

2.5089 

7.8375 

9.00 

1.2040 

1.8429 

2.2303 

2.5089 

8.8375 

10.00 

1.2040 

1.8429 

2.2303 

2.5090 

9.8375 

OB 

1.2040 

1.8429 

2.2303 

2.5090 

00 

P*  » 

.90 

fV 

2 

3 

4 

5 

CD 

.25 

.1219 

.1410 

.1440 

.1445 

.1446 

.50 

.3359 

.3785 

.3892 

.3927 

.3946 

.75 

.5399 

.6085 

.6289 

.6372 

.6446 

1.00 

.7314 

.8306 

.8624 

.8768 

.8946 

1.25 

.9077 

1.0437 

1.0892 

1.1109 

1.1446 

1.50 

1.0656 

1.2457 

1.3079 

1.3385 

1.3946 

1.75 

1.2023 

1.4342 

1.5169 

1.5584 

1.6446 

2.00 

1.3158 

1.6064 

1.7140 

1.7689 

1.8946 

2.25 

1.4053 

1.7594 

1.8966 

1.9679 

2.1446 

2.50 

1.4724 

1.8905 
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